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Abstract: In models with partially composite quarks, the couplings of quarks to the Z 
boson generically receive non-universal corrections that are not only constrained by elec- 
troweak precision tests but also lead to flavour-changing neutral currents at tree level. The 
impact of these flavour-changing couplings on rare K and B decays is studied in two-site 
models for three scenarios: an anarchic strong sector with two different choices of fermion 
representations both leading to a custodial protection of the Z — >• bb coupling, and for a 
strong sector invariant under a U (2) 3 flavour symmetry. In the complete numerical analy- 
sis, all relevant constraints from AF = 2 processes are taken into account. In all scenarios, 
visible effects in rare K and B decays like K —> ttvP, B s — > fJ, + fJ,~ and B — > K* fi + /j,~ are 
possible that can be scrutinized experimentally in the near future. Characteristic corre- 
lations between observables allow to distinguish the different cases. To sample the large 
parameter space of the anarchic models, a new method is presented that allows larger 
statistics than conventional approaches. 



1 Introduction 



The discovery of a Higgs boson with Standard-Model-like properties faces us with a severe 
fine-tuning problem. A straightforward solution to this puzzle would be a composite nature 
of the Higgs scalar. To construct viable models with a composite Higgs boson, fermion 
masses have to be generated without at the same time violating the stringent bounds on 
flavour-changing neutral currents (FCNCs). A compelling mechanism to achieve this goal is 
partial compositeness, where the SM fermions are assumed to couple linearly to heavy, com- 
posite fermions with the same quantum numbers [1]. Fermion masses are generated at low 
energies in a see-saw-like manner and the mass eigenstates are superpositions of elementary 
and composite fields, the degree of compositeness being larger for heavier fermions. This 
mechanism, explaining the fermion mass hierarchies by hierarchical composite-elementary 
mixings, in turn leads to a suppression of FCNCs even if the strong sector is completely 
flavour- anarchic [2-4]. Still, barring cancellations, the strong bounds from CP violation 
in kaon mixing require the resonances of the strong sector to be in excess of 10 TeV, in 
tension with naturalness [5, 6]. Consequently, various mechanisms have been suggested to 
further suppress flavour violation, one possibility being a flavour symmetry under which 
the strong sector is invariant, only broken by the composite-elementary mixings [7-1 1] 1 . 

Another challenge of these models are electroweak precision tests. To avoid excessive 
contributions to the T parameter, the strong sector should be invariant under the custodial 
symmetry SU(2)l x SU(2)r x U(1)x, containing SM hypercharge as Y = T^r + X [16]. 
Moreover, the couplings of the Z boson to fermions, measured with high precision at LEP, 
are modified by two effects. First, since the heavy fermions come in vector-like pairs, 
after electroweak symmetry breaking (EWSB) the SM fermions mix with states that have 
different electroweak quantum numbers. Second, since the composite states are charged 
under the enlarged custodial symmetry, the mixing of vector resonances associated to 
SU(2)l and SU(2)r leads to additional corretions after EWSB. Both effects can also lead 
to flavour-changing Z couplings after going to the mass eigenstate basis. The resulting tree- 
level FCNCs put additional constraints on these models. The strongest flavour-conserving 
constraint, which comes from the Z coupling to left-handed b quarks, can be avoided by 
making the strong sector (approximately) invariant under a discrete symmetry, which can 
be achieved by judicious choice of fermion representations [17]. This choice will affect 
flavour-changing Z couplings as well, leaving a characteristic pattern of effects than can be 
probed by FCNC measurements. 

The experimental prospects for probing these flavour-changing Z couplings in the near 
future are excellent. In the b — > s sector, the recent measurement of the decay B s — > // + fj,~ 
[18] has excluded sizable scalar current contributions to this mode and has started to 
become sensitive to Z-mediated effects. At the same time, the ongoing measurements 
of angular observables and CP asymmetries in the B — > K* fi + [i~ decay will allow to 
disentangle chiralities and phases of the flavour- changing Z couplings. In the s — > d sector, 
improved experiments probing the decays Kl — > iftvv and K + — > tt + vv are imminent. 

Alternative solutions include flavour symmetries broken also in the strong sector [12-14] or an extension 
of the (flavour-blind) global symmetry of the strong sector [15]. 
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The goal of this paper is a detailed numerical study of the effects of flavour-changing Z 
couplings in models with partial compositeness. The focus will be on two different choices 
for the heavy fermion representations, both designed to avoid excessive contributions to 
the Zbb coupling, as well as on two different flavour structures: first, a scenario with an 
anarchic strong sector; second, a scenario with a strong sector invariant under a C/(2) 3 
flavour symmetry [10, 19, 20], minimally broken by the composite-elementary mixings. 
The basic structure of these models and the order of magnitude of the relevant effects have 
recently been analyzed in ref. [20]. Here, they will be scrutinized by a complete numerical 
study. The main novelties of this paper are 

• a first full numerical analysis of the custodially protected flavour-anarchic model 
where the left-handed up- and down-type quarks couple to two different composite 
fermions (to be called bidoublet model below), 

• a full numerical analysis of the custodially protected flavour- anarchic model where 
the left-handed quarks couple to a single composite fermion (to be called triplet model 
below), which has already been analyzed numerically in a similar framework in the 
literature [21, 22], using updated flavour constraints (e.g. the B s mixing phase), 

• a novel method to sample the parameter space of partially composite models with 
flavour anarchy that allows much larger statistics than conventional methods, 

• a first full numerical analysis of the model with a f/(2) 3 flavour symmetry in the 
strong sector [10, 20], minimally broken by the composite-elementary mixings of 
right-handed quarks. 

The paper is organized as follows. In section 2, the definition of the framework and the 
different scenarios will be recalled. In section 3, analytical expressions for the Z couplings in 
the different scenarios will be derived to obtain a qualitative understanding of the effects 
in the numerical analysis. Section 4 will briefly review the relevant FCNC observables. 
Section 5 is devoted to the numerical analysis. The summary and conclusions will be 
presented in section 6. The novel method to scan the parameter space of anarchic models 
is discussed in appendix A. 

2 Setup 

A simple effective framework to study the phenomenology of the models of interest is given 
by the two-site approach of ref. [23]. In this framework, one considers one set of fermion 
resonances with heavy Dirac masses as well as spin-1 resonances associated to the global 
symmetry SU (3) c x SU(2) l X SU(2)r x U (1)x ("heavy gauge bosons"). This approach can 
be viewed as a truncation of 5D warped models, taking into account only the lightest set of 
KK states. This approximation is particularly justified in the case of tree-level amplitudes, 
which are the only ones relevant for this study. As a further simplification, all vector 
resonances will be assumed to have a common mass and coupling, m p and g p , respectively. 
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As mentioned in the introduction, we will consider two different choices for the heavy 
fermion representations, dubbed triplet and bidoublet models. Starting with the triplet 
model 2 , its fermion Lagrangian reads 



C m = -tilVm^V] - tx^m^W] - tr[R n m%R ,j ] 
Y i hT\L i L UR^ + F ij tr[ft V L R^\ + h.c. 



(2.1) 
(2.2) 



where capital letters denote composite fields, which transform under SU(2)l x SU(2)r x 
U(l) x as 



L = (Q Qf) 
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While the strong sector C m + C Y is invariant under the full custodial group, the composite- 
elementary mixings break SU(2)r x U(1)x explicitly. 
The Lagrangian of the bidoublet model 3 reads 



C m = -tr^mgX] " U ir 4ui + (U D) 
YHtx[L\jU\ L U 3 R + h.c.) + (U -+ D) 



^ mix = + 4Jji*k + (U,u^ D, d) 



where 



Lu = (Qu Q'u) 
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(2.4) 
(2.5) 
(2.6) 

(2.7) 
(2.8) 



An important new feature with respect to the triplet model is that the left-handed up- 
and down- type quarks now couple to two different composite fields, so that the degrees 
of compositeness of the left-handed top and bottom quarks, for instance, are unrelated. 
This is crucial not only for the suppression of the corrections to the Z — > bb coupling, as 
discussed in section 3, but is also necessary in models with a flavour symmetry only broken 
by the left-handed composite-elementary mixings [9]. 

For the flavour structure, we consider two possibilities: anarchy and C/(2) 3 . In the 
anarchic case, the strong sector has no flavour symmetry and the quark mass and mix- 
ing hierarchies are generated by hierarchies in the composite-elementary mixings. In the 



2 The triplet model was called TS10 (the two-site model with the fermions fitting into a 10 of SO(5) x 
U(l)x) in ref. [24]. A variation of the triplet model is obtained by having the Ur mix with an additional 
(1; 1)2/3 instead of the U from the triplet (see e.g. [21, 22, 25, 26]). While the protection of Z couplings is 
the same as in the triplet model, due to the larger number of parameters, FCNCs effects tend to be even 
bigger. We will stick to the minimal case without the singlet in the following. 

3 The bidoublet model was denoted TS5 in refs. [24, 27] (the fermions fit into two 5 of SO (5) x U(l)x). 
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bidoublet model, the Yukawa couplings Yff D are anarchic matrices, while the composite- 
elementary mixings A^ can be chosen to be real and diagonal. In the triplet model, a 
basis where A^- are real and diagonal can be chosen as well, although SU(2)r is no longer 
manifest in this basis and one effectively gets different Yukawas in the up- and down-type 
sectors, related by Yy = Y^U^, where U is a unitary matrix. For simplicity, the Dirac 
mass terms will be assumed to be proportional to the identity in the following 4 . 

In the U (2) 3 case, the strong sector is assumed to be invariant under a flavour symmetry 
under which the first two generation fields transform as doublets and the third generation 
as singlets [19]. The symmetry is broken by the left-handed or right-handed composite- 
elementary mixings [10, 20]. The breaking is assumed to be minimal, i.e. using the smallest 
set of small symmetry breaking spurions that is required to generate the observed quark 
masses and mixings. The case where the breaking occurs in the right-handed mixings is 
called left-handed compositeness. In this case, the Dirac mass matrices, Yukawa matrices 
and left-handed mixings are all of the form diag(a, a, b) with real entries. The form of 
the right-handed mixings compatible with the minimal breaking of C/(2) 3 can be found in 
the appendix of ref. [20]. Analogously, in right-handed compositeness, the right-handed 
mixings are diagonal while the left-handed ones contain the U(2) 3 breaking spurions. 

A further attractive possibility is that the strong sector is fully flavour- invariant, i.e. 
possesses a U(3) 3 symmetry [7-9]. In this case however, flavour-changing Z couplings do 
not arise at tree-level at all, so it will not be considered in the following. 

3 Z couplings 

In ref. [17] it has been shown that two discrete subgroups of the custodial group can 
protect the Zbb coupling from excessive tree-level corrections, as required by electroweak 
precision measurements: the Plr symmetry, which exchanges SU(2)l and SU(2)r, and 
the Pc symmetry, under which T^l r —T3L,r- In the triplet model, the d L mix with 
the Plr eigenstates B l , while the u R mix with the Pc eigenstates U % . Consequently, 
the (flavour-conserving and flavour-changing) Z couplings to left-handed down-type and 
right-handed up-type quarks at zero momentum receive no correction at tree-level. In the 
bidoublet model, the right-handed quarks mix with the Pc and Plr eigenstates U and 
D, so the right-handed Z couplings are not corrected. In the case of left-handed mixings, 
the situation is slightly more involved compared to the triplet model, since the left-handed 
elementary doublet now mixes with two different fields, Q u and Qd, via the mixing terms 
Xlu and XlcL- While the T' and B (in the notation of eq. (2.8)) are Plr eigenstates, the 
T and B' are not. Consequently, the Zd L d J L couplings receive corrections involving XlcL-, 
but not Xlu- This is enough to protect the Zbi^L coupling, since a hierarchy \Ld *C Xlu 
is natural in view of the lightness of the b quark compared to the top quark. 

To obtain approximate analytical expressions for the modified Z couplings, let us start 

4 Relaxing this assumption will not lead to qualitatively new effects, but might allow larger effects in 
view of the enlarged parameter space. 
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Table 1. Tree-level contributions to the modified Z couplings as defined in eq. (3.1) in the triplet 
and bidoublet models. The Ax are defined in eq. (3.3). 

with a simplified one-generation model. Writing the Z couplings to quarks as 

-q^ [(T 3L ~ Qs 2 w + 5g Lq )P L + {-Qs 2 w + 5g Rq )P R ] q Z^ , (3.1) 



g_ 



the contributions due to fermion and gauge boson mixing can then be written as in table 1, 
where 

sx = Xx/^m 2 x + X 2 x (3.2) 

is the degree of compositeness and 

A - v2Y v A - v2Y » A - v2Y ° A - f3 3) 

Alu "2^' ALd "2^' AM "2^i' As "4^f (3 - 3) 

In the three-generation case, (3.3) remains valid for the flavour-diagonal couplings 5 , in 
particular for the Z — >■ bb coupling, while the flavour-violating couplings are dependent 
on the flavour structure. In anarchy, all the flavour-changing couplings allowed by the 
Plr and Pq symmetries are generated, while in C/(2) 3 , the right-handed flavour-changing 
couplings are strongly suppressed. In fact, in t7(2) 3 with right-handed compositeness, no 
flavour-changing couplings are generated at all at tree- level [20]. A generic prediction of 
minimally broken t7(2) 3 is that the amplitudes of 12-transitions are aligned in phase with 
the SM [19]. 

With these considerations, one can already recognize a pattern among the flavour- 
changing Z couplings allowed by the different electroweak and flavour structures that can 
potentially serve to distinguish them experimentally. Table 2 lists the presence or absence 
of the flavour-changing couplings relevant for K, Bd tS and D decays for the models discussed 
so far. 

One can now go one step further and estimate the allowed size of the flavour-changing 
couplings in the various models, to understand the size of the effects encountered in the 
full-fledged numerical analysis of section 5. To this end, it turns out to be convenient to 
consider rescaled couplings 

€n= 8 $f L ^l d ,Rd^ (3-4) 

where £jj = VtjV t *. This definition is chosen so that 0(1) values of the 5 l [ R correspond 
roughly to the size of the loop-induced flavour-changing Z coupling in the SM (and thus can 
lead to visible effects in FCNC observables, see section 4) and for real 8£ R , the couplings 
are aligned with the SM contribution. 



In anarchy, this is true if one interprets the Yu,d as "average" anarchic Yukawas. 
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Table 2. Pattern of effects in flavour-changing Z couplings relevant for K (s — > d), B (b — > s, d) 
and -D (c — » u) physics. R denotes contributions aligned in phase with the SM, (D contributions 
with a new CPV phase and empty cells refer to vanishing or strongly suppressed contributions. 

Anarchic triplet model Here the flavour-changing couplings can be estimated by re- 
placing s 2 R — > s^s^ in table 1. Using the well-known approximate relations 6 between the 
degrees of compositeness, the CKM angles and the quark masses, 



Y* s Ld s Rd i 



Vij ~ s l Ld /s 3 Ld , 



(3.5) 



and defining xt = SLt/sRt ~ 0(1) > one can write the right-handed flavour-changing cou- 
plings to down-type quarks as 



Sg 



Vd^VdJ 
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\A Rd + A g ) 
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(3.6) 



(3.7) 



Consequently, one can expect sizable effects in K decays, while effects in B or B s decays 
should be smaller. Note that the size of these effects is not constrained by Z — > bb, which 
mostly constrains the left-handed Zbb coupling (cf. [20, 28]). 

Anarchic bidoublet model Replacing accordingly s\ — > s^s^ in table 1, one can es- 
timate the left-handed flavour-changing couplings in terms of the coupling to b quarks 

as 



33 

Ld ■ 



Since 5gf^ d is constrained by Z — )■ bb to be less than about 5 x 10 4 , one finds 



I s-12,13,23 1 n -i 
\0 L | ~ 0.1 



5 x 10 



-4 



(3.8) 



(3.9) 



Consequently, one expects effects in K physics to be smaller than in the triplet model, but 
could have sizable effects in B physics. 



Replacing all functions of Yukawa couplings by appropriate powers of an "average" Yukawa Y* 
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U(2) 3 bidoublet model In this case, the left-handed flavour-changing couplings to 
down-type quarks can be estimated as 



33 



33 



(3.10) 



where is a complex parameter with magnitude of O(l). Similarly to the anarchic bidou- 
blet model above, one thus obtains 



.-12,13,23 1 



0.1{|r b | 2 ,|r b |,|r b |} 



5g 33 



Ld 



5 x irr 4 



(3.11) 



For \rt>\ ~ 1, one could thus have visible effects in both K and B decays, while for \r^\ < 1 
(> 1) effects in B (K) physics should be larger. 



4 Probes of flavour-changing Z couplings 

The flavour-changing Z couplings contribute to FCNC processes mediated by Z exchange in 
the SM, i.e. rare K, B, B s and D decays. Below some formulae for the relevant observables 
are collected. 

4.1 s d 

The branching ratios of the charged and neutral K — > ttvv decays can be written as 

2 / r> ti- v w 2" 



BR(K + tt + uu) = k+(1 + A 



EM J 



— a* — ) +^V) + — # — ) 



BR(K L 



IT VV) 



Am(ei2^) V 



(4.1) 
(4.2) 



where X K = X SM + (5 L 2 + 5^) with X SM = 1-469 ± 0.017 [29]. The quantities «+, k l and 
P(u,c) can be found e.g. in [29]. The experimental measurement and SM predicitions read 
[29]' 



7T VV 



cxp 



;i7.3i£jo x 10 



-n 



br(k a 

BR(K + -> tt + vv) sm = (7.81 ± 0.80) x 10~ n 
BR{K L -> 7t ^)sm = (2.43 ± 0.39) x 10~ n 



(4.3) 
(4.4) 
(4.5) 



while the experimental bound on Kl — > tt°vv has not yet reached the model-independent 
theoretical upper bound [30]. 

Another relevant process is the decay Kl — > H + H~ ■ Unfortunately, it is plagued by 
long-distance contributions which are hard to calculate. The short-distance part can be 
written as 



BR{K L 



M A* JSD 



A 5 



(4.6) 
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^23 _ x23 2 
1+ *L *R_ 



(4.8) 



where Y K = Y SM + {5f - 5^) with Ism = 0.98 ± 0.02 and the quantity P c can be found 
in [31]. An estimate of the long-distance contributions allows to extract an experimental 
upper bound on the short-distance part [32] 

BR(K L -> ^ + /Osd < 2 - 5 x 10_9 • ( 4 - 7 ) 

Since Kl — > depends on the difference of left- and right-handed couplings, while the 

K — > nvi> decays depend on their sum, the constraints are complementary. 

4.2 b -> s 

The b — > s transition is the one best constrained experimentally due to a multitude of 
inclusive and exclusive decays sensitive to it. For our purposes, the most relevant decays 
are B s — > [i + fi~ and B — > K*fi + ^~. 

In the absence of scalar currents, the branching ratio of the B s — > /i + \i~ decay can be 
written as 

BR(B,-»/xV) 
BR{B S ^ /i+/i-) S M 

where Ism was given after eq. (4.6). For the flavour-averaged branching ratio, one currently 
has [18, 33] 

BR(-Bs -> /x + /O CX p = (3.5±i:|) x 10- 9 , (4.9) 
BR(B S -> /x + /Osm = (3.23 ± 0.27) x 1(T 9 . (4.10) 

Eq. (4.8) is also valid for the decay B^ — > H + H~ with the replacement 5^ ~~ ^ &x • 

The B — > K*fi + n~ decay allows access to many observables sensitive to new physics 
by means of an angular analysis. Since only modified Z couplings are of interest here, 
only the region of large dilepton invariant mass, q 2 > 14.18 GeV 2 , will be considered. This 
region has two advantages. First, magnetic dipole operators, which can also be generated 
in the models at hand but are not taken into account here, do not contribute in this region. 
Second, this region is theoretically cleaner since non-factorizable effects are under better 
control [34] and since the form factors can be calculated by lattice QCD. 

While the precise expressions for the B — > K*fi + fi~ observables can be found e.g. in 
ref. [35] , below approximate numerical expressions for the central values of theory param- 
eters are reported to understand the sensitivity of the observables to the Z couplings (in 
the numerical analysis the exact formulae are used). For the branching ratio, one finds 



1 + 1.3 Re(<5f)- 0.91 Re(<5£ 



+ 0.63 (|<5i 3 | 2 + |4 3 | 2 ) " 0.91Re(5f4 3 *) • (4.11) 

Among the angular observables, the forward-backward asymmetry Ap^, the angular CP 
asymmetry Ag and the CP-averaged angular observable 5*3 are particularly interesting, 
since they can be extracted from a one-dimensional angular distribution. Expanding their 
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precise expressions to first order in the couplings S^ R , one obtains for them 



A FB (B ET>+^-) SM 



1 - 0.12Re(<5f ) - 0.93 ReO^ ) , (4.12) 



l-1.4Re(<$£), (4.13) 



S 3 (# -> K* M +/i-) SM 

A)(£ -»• K*fi + fi-) « -0.31 Im(Sf ) , (4.14) 

valid for the high g 2 region, g 2 > 14.18 GeV 2 . Importantly, the observables S3 and are 
only sensitive to the right-handed coupling. Experimental averages and SM predictions 
for all the observables can be found e.g. in ref. [36]. 

While not accessible at the LHC, future B factories will look for the decays B — > K^*'uv 
(and possibly also the inclusive decay). Being free from photon penguin contributions, also 
these decays are clean probes of modified Z couplings. Defining 



IW + W -Re (5?5f*) , 

v= ,^,9 , ( 4 - 15 ) 



A SM ' \Sf\ 2 + \S%\ 2 ' 

where Xsm was given in section 4.1, the branching ratios and the angular observable Fl 
in B — > K*vv can be written as [37] 

BR(5 -> Kvv) , . , BR(5 -»■ iT*i/p) , . , 

BBiB^X**) _ ll + 0M yi - /' +2 "' ■ 



BR(B -> X,Ki/)sM ft,(B -> K*kP) S m (1 + 1.31 ij) 

A combined analysis of these observables would thus allow to disentangle left- and right 
handed couplings. 



5 Numerical analysis 

After estimating the size of the flavour-changing couplings in section 3 and defining the 
relevant experimental observables in section 4, we are now ready for a numerical analysis of 
the effects in the various models considered. Such analysis is particularly challenging in the 
anarchic case in view of the large number of free parameters and the necessity to reproduce 
the known quark masses and mixings. To this end, a novel approach was used in this study, 
based on a Markov chain Monte Carlo and described in more detail in appendix A. Here it 
suffices to say that the result corresponds to a sample of the parameter space with certain 
ranges allowed for the free parameters. Throughout the numerical analysis, a common 
mass and coupling, m p and g p , was assumed for the vector resonances. In the anarchic 
case, as mentioned above, the further simplification of assuming the heavy fermion Dirac 
mass matrices to be proportional to the identity was made. The parameters were then 
allowed to float in the ranges 

9p e[3,5], y« D €[0,47r], arg^) G [0, 2vr] , (5.1) 

m p G [3, 5] TeV , rrty G [0.5, 3] TeV , (5.2) 
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Figure 1. Right-handed flavour-changing Z couplings in the anarchic triplet model relevant 
for s s — > d (left) and b — > s transitions (right), normalized as in eq. (3.4). The gray regions show 
the combined experimental constraints at 95% C.L. In the left-hand plot, the dashed circle shows 
the constraint from BR(X + — > n + vv) alone. In the right-hand plot, the dashed ring shows the 
constraint from BR(B S — > h + /j,~) alone. 

where stands for any heavy fermion mass, with all the additional parameters uncon- 
strained. Since the goal of the present analysis is to obtain the maximum allowed deviations 
of the flavour- changing Z couplings from the SM, the parameter ranges were chosen to allow 
for sizable effects. The lower bound on m p was chosen to fulfill the constraint from the S 
parameter, while the upper bound on g p was chosen to always have / = rn p /g p > 600 GeV. 
Note that, contrary to ref. [20], the Ytj.d and the composite fermion masses where treated 
as unrelated. 

Throughout the analysis, the relevant constraints from AF = 2 processes, namely 
the mass differences and mixing phases of the Bj and B s systems and the CP violating 
parameter |e^|, have been required to be within their experimentally allowed bounds (for 
details of the procedure see appendix A). As is well known, the |ex| constraint is particularly 
hard to fulfill in the flavour anarchic case. The approach of this study is to accept the 
implied fine-tuning in the parameter space and to study the maximum allowed size of 
flavour-changing Z couplings. In the U(2) 3 case, this problem is absent. 

5.1 Flavour-changing Z couplings 

We start by looking at the size of s — >■ d and b — > s couplings in the various models. 
Figure 1 shows <5^ 2 and 5jf in the anarchic triplet model, compared to the experimentally 
allowed ranges. All points satisfy the constraints from AF = 2 observables and Z — > bb. In 
the left-hand plot, the gray region shows the combined constraint from BR(i<C + — > n + vv) 
and BR(Kl — > /i + /i~)sD (treating eq. (4.7) as a 68% C.L. bound), while the dashed circle 
shows the constraint from BR(ivT + — > ir + vv) alone. In the right-hand plot, the gray region 
shows the combined constraint from all relevant B — > K*fi + ^i~ observables at high q 2 and 
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Figure 2. Left-handed flavour-changing Z couplings in the anarchic bidoublet model relevant 
for s — >■ d (left) and b — > s transitions (right), normalized as in eq. (3.4). 

from BR(i? s — > while the dashed ring shows the constraint from BR(1? S — > /z + /i~) 

alone 7 . As expected from the estimates in section 3, the effects in K physics are sizable, 
while the effects in B physics are smaller (but still sizable compared to the more stringent 
experimental bounds). Note that the tilted main branch in the left-hand plot corresponds 
to almost real 5g^ d , which is preferred due to the interplay between AS = 1 and AS" = 2 
amplitudes and the strong constraint on \ ex\ ■ The tilt is caused by dividing out the complex 
SM CKM factors in eq. (3.4). 

Figure 2 shows analogous plots in the anarchic bidoublet model, this time showing 
the left-handed couplings. Also here, the estimates of section 3 are confirmed: while the 
effects in the s — > d coupling are smaller than in the triplet case, the effects in b — > s can 
be sizable. In fact, both B —> K*fi + fj>~ and B s -4 fi + already cut into the model's 
parameter space. 

Figure 3 shows the size of the left-handed b — > s couplings in the bidoublet model with 
U (2) 3 (left-handed compositeness) . Just as in the anarchic case, one can obtain large effects 
in B physics with a new CP violating phase. The main difference is that the size of the 
imaginary part is now limited. The reason is that, for the parameter r\, in (3.10) to deviate 
from 1 (in which case the Z coupling would be aligned in phase with the SM) requires a 
relatively large degree of compositeness for the first two generations of left-handed up-type 
quarks, leading to a sizable correction to the hadronic width (see the discussion in the 
appendix of ref. [20]). The effects in the s — > d couplings are similar in magnitude to the 
anarchic case, but now Im^J 2 ) = 0, which is why the corresponding plot is omitted. 

7 For details on how these constraints are obtained numerically, the reader is referred to ref. [36]. 
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Figure 3. Left-handed flavour-changing Z couplings in the bidoublet model with U(2) s left- 
handed compositeness relevant for b — > s transitions, normalized as in eq. (3.4). 



5.2 Predictions for rare decays 

Having established the size of the modified Z couplings, we can now proceed to analyze 
their impact on the FCNC observables 8 . 

Figure 4 shows the correlations between BR(K + — > i: + i>v) and BR(Kl — > it vis), 
between BR(iT + — >• ir + vv) and RR{Kl — > /i + ^~)gD, between BR(5 S — > /i + ^~) and 
BR(B d -4 between BR(B S -4 and BR(B -4 K*fi + fi-) (in the high q 2 re- 

gion), between the B -4 K*fj, + fi~ observables £3 and Ag as well as between BR(B —> K*vi>) 
and BR(B — > Kvu) in the anarchic triplet model, the anarchic bidoublet model and the 
bidoublet model with U(2) 3 left-handed compositeness. In all plots, the gray regions cor- 
respond to the experimental la ranges. In the case of BR(Kl — > H + H~)st>, the horizontal 
gray line shows the bound in eq. (4.7). In the upper left plot, the diagonal line shows the 
model-independent Grossman-Nir bound [30]. The black cross in each plot shows the SM 
central values and uncertainties. All the new physics points are obtained for central values 
of the theory parameters, so are understood to carry at least the relative uncertainty of 
the SM point. 

In the case of the anarchic triplet model, represented by the yellow points in figure 4, 
the following observations can be made. 

• K + —> ■k + vv and Kl — > already cut into the model's parameter space and 
BR(Kl — > iflvv) can be enhanced to five times its SM value. 

• BR(K + —7- TT + vi>) and BR(Kl — > /U + /i~)sD are anticorrelated. This is because only 
the right-handed flavour-changing Z coupling contributes. 

8 Additional contributions to AF — 1 amplitudes from the exchange of neutral vector resonances are 
subleading and thus neglected. For AF = 2 constraints, they are instead important and taken into account. 
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Figure 4. Correlations between observables in rare K and B decays in the anarchic triplet model 
(yellow/light gray), the anarchic bidoublet model (red/medium gray) and the bidoublet model with 
[/(2) 3 left-handed compositeness (blue/dark gray). The shaded regions indicate experimental la 
ranges. The black cross shows the SM predictions with theory uncertainties. 
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• The bound on Kl — > fj> + /i excludes a suppression of BK(K + — > ir + vv) by more 
than 50%. 

• Both B s — > fi + fj>~ and Bd — > h + /j>~ can be suppressed or enhanced by up to 50%. 
The precision of the B s — > (j7 ' \i~ measurement is now entering the region accessible 
in the model. 

• BR(f? s — > and BR(£? — > K*n + fi~) are perfectly correlated. Prom equations 
(4.8) and (4.11) one can see that such correlation is present if only left-handed or 
only right-handed couplings (as in this case) contribute, while deviations are possible 
if both left- and right-handed couplings are present. 

• The CP asymmetry Ag in B — > K*fi + ^L~ at high q 2 can reach at most 5%; the CP- 
averaged observable S3 can deviate by ±10% from its SM expectation. At present, 
the theory uncertainty due to form factors is too large to constrain the model. 

• BR(S — > Kuv) and BR(5 — > K*uv) are anticorrelated, again because only the 
right-handed flavour-changing Z coupling contributes (cf. [37-40]). 

Note that although the experimental measurement of the B — > K*fi + fi~ branching ratio 
is already quite precise, the large SM uncertainty limits its constraining power at present. 
The situation will change once high-(7 2 form factors will be computed on the lattice. 

The red points in figure 4 show the analogous correlations in the anarchic bidoublet 
model. One can make the following observations. 

• The effects in K + — > tt + uv and Kl — > t^vv are almost unconstrained by experiment 
as of yet, but enhancements by factors of 2 are possible. 

• BR(A' + — > tt + uv) and T5H(Kl — > n + fi~)sD are now perfectly correlated rather than 
anticorrelated since only the left-handed flavour-changing Z coupling contributes. 

• The measurement of B s — > already constrains the model signifcantly, while 
effects in B^ — > fi + fi~ are at most ±50%, not yet constrained by experiment. The 
two decays are uncorrected. 

• BR(£? S —J- [i + vs. BR(i? — > K* n + fj>~) shows almost the same correlation as in the 
triplet model. This is because now, only left-handed couplings contribute. 

• The observables S3 and Ag in B — > K* fi + ^ are SM-like since there are no right- 
handed flavour-changing Z couplings. 

• Also BR(i? —7- Kvv) and BR(i? — > K*vv) are now correlated rather than anticorre- 
lated, because the effects are due to left-handed couplings. 

Finally, the dark blue points in figure 4 show the correlations for the U (2) 3 bidoublet 
model with left-handed compositeness. There are two main differences with respect to the 
anarchic bidoublet model. 



- 15 - 



• K + —7- 7r + zy^ and Kl — > irw are perfectly correlated since the amplitude is aligned 
in phase with the SM. 

• Also B s — > fM + fj," and — > fj> + ^ are perfectly correlated since b — > d and b — > s 
amplitudes are universal, up to the same CKM factors as in the SM. 

Both these correlations are universal predictions of models with Minimal Flavour Violation 
(MFV) [41, 42] and of models with a minimally broken U(2) 3 symmetry (as in this case). 

6 Conclusions and Outlook 

In models with partial compositeness, flavour-changing Z couplings can lead to visible 
new physics effects in rare K and B decays. The chirality, phase and size of these flavour- 
changing couplings depend on the one hand on the mechanism imposed to keep the flavour- 
conserving Z bb coupling under control and on the other hand on the presence or absence 
of a flavour symmetry imposed to ameliorate the €k problem. In this paper, these effects 
have been analyzed numerically in three scenarios, two flavour-anarchic ones with different 
fermion content (dubbed triplet and bidoublet models) and a model with a U(2) 3 flavour 
symmetry, minimally broken only by the composite-elementary mixings of right-handed 
quarks. In all cases, a two-site Lagrangian was used with some simplifying assumptions, 
e.g. a common mass and coupling for all the vector resonances. A new method was used 
to sample the highly-dimensional parameter space of these models. 

The main result of this study is figure 4. It shows the size of and correlations between 
the most relevant rare K and B decays. The observables that have been chosen are not 
sensitive to magnetic dipole operators that are also generated in these models, so they can 
be considered clean probes of flavour-changing Z couplings. 

Interestingly, in all the scenarios, all the branching ratios considered can exhibit visible 
effects at current and planned experiments such as LHCb, Belle-II, KOTO and NA62. 
Particularly large effects, already constrained by existing experimental results, are possible 
in K decays in the anarchic triplet model and in B decays in the case of the bidoublet 
models. Clear-cut correlations between various observables can be used to distinguish the 
three models on the basis of the pattern of observed deviations from the SM expectations. 

This study could be extended in various ways. Going beyond the simple two-site 
picture and adopting more complete frameworks, e.g. in a 5D picture, additional relations 
between the parameters could modify the maximum size of the effects found here. In 
addition, a number of potentially relevant constraints has not been taken into account in 
this analysis, namely ones related to loop-induced observables. The rationale is that loop 
effects are more dependent on the details of the theory, e.g. the localization of the Higgs 
in the extra dimension in a 5D picture (see e.g. [26, 43, 44]). In concrete frameworks 
however, such effects can be calculable. The most important loop-induced quantities are 
the T parameter and magnetic dipole operators contributing to FCNCs and electric dipole 
moments. Their study is left to a future publication. 
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A Details on the numerical analysis 
A.l Scanning procedure 

The numerical analysis of partial compositeness models with anarchic Yukawa couplings is 
challenging since the quark masses and mixings depend on a large number of free parame- 
ters, making a brute-force scan of the highly-dimensional parameter space unfeasible. The 
conventional approach to tackle this problem is to exploit approximate analytical relations 
between the Yukawas, degrees of compositeness and the masses and mixings [21, 45, 46]. 
This approach has several drawbacks. First, it is difficult to obtain the correct CKM phase, 
so that a partially brute-force approach is necessary. Second, the approximate analytical 
expressions become very complicated if the left-handed quark doublet couples to two dif- 
ferent fields, as is the case in the bidoublet model. Third, once a set of viable parameter 
points is obtained, imposing additional constraints, like Z — > bb or AF = 2 observables 
(especially \€k\), one looses a large fractions of the points. 

To solve these problems, this study uses an alternative approach based on Markov 
Chain Monte Carlo (MCMC) with the Metropolis-Hastings algorithm. This method is a 
very simple tool to sample a highly-dimensional probability distribution, given a likelihood 
function and a prior distribution for the input parameters. While it is frequently used 
for Bayesian parameter estimation, here a different property is exploited: the algorithm 
automatically samples the region in parameter space where the likelihood is high, but 
(in the stationary limit, i.e. with infinitely long chains) is supposed to sample the entire 
relevant parameter space. 

The procedure is as follows. One defines the likelihood as L = e - * 2 ^/ 2 , where the x 2 
function depends on all the model input parameters 8 (including the real and imaginary 
parts of the Yukawa matrices, the composite-elementary mixings, etc.) and contains the 
experimentally measured quark masses and CKM angles. One now starts a chain with step 
size chosen to obtain an acceptance rate around 20%. That is, after a burn-in period, out 
of 5 generated points, 1 point is viable. The chain has to be long enough to obtain an 
acceptable coverage of the parameter space. Several chains with random initial values can 
be combined to make sure disjoint minima are not missed. 

One can now go one step further and include additional constraints in the x 2 ■ This 
was done in the numerical analysis of section 5 for the Z — >■ bb partial width and the Z 
hadronic width as well as for the AF = 2 observables AM^, AM S / AM^, S^k S i 0s an d £k- 
Apart from drastically reducing the number of points to be generated - since constraints are 
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Table 3. Experimental values and hadronic parameters used for the AF — 2 constraints. 

automatically fulfilled - this has an additional advantage: in the conventional approach, the 
CKM angles and phase to be used in the first step are always the same, while the AF = 2 
processes calculated in the second step actually change the preferred CKM parameters. If 
AF = 2 constraints are included alongside CKM angles obtained from tree-level processes, 
the impact of new physics effects in loops on the CKM determination is automatically 
taken into account. For example, a sizable constructive new physics contribution to the 
Bd mixing phase leads to a lower value of p in the determined point, so that S^k s agrees 
with the experimental value. 

While the acceptance rate can be kept high by adjusting the step size, the downside 
of including additional constraints is an increase in the required burn-in time as well as 
a reduced coverage of the parameter space for a given chain length. Consequently, one 
needs a larger number of points to obtain a representative sample. For each of the plots in 
section 5, approximately 10 6 viable points, consisting of 5-10 chains, were used. 

It should be stressed again that the only purpose of the MCMC approach is to sample 
the regions in parameter space which have the right quark masses and mixings and which 
fulfill the flavour cosntraints. The purpose is not to obtain posterior probability distri- 
butions for the input parameters. This would be a futile task since the number of input 
parameters far exceeds the number of constraints. The dependence on the choice of priors 
and the limited coverage are thus much less problematic than they would be in the case of 
a probabilistic analysis. 

A. 2 Input parameters and constraints 

For the inclusion of the AF = 2 constraints, the Wilson coefficients of the operators 

Qv L = (47 M 4)K7 M 4) , Qv R = (47^4X47^4) , (A-l) 

Qv R = (47"4)(47 M 4) , Qs R = (44X44) . (a.2) 

with ij = 21, 31, 32 are evaluated at tree level, including contributions from neutral coloured 
and uncoloured vector resonances. The renormalization group evolution from the scale 
m p to the appropriate low-energy scale follows ref. [57]. The input parameters used are 
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collected in table 3. The CKM Wolfenstein parameters correspond to a fit to tree- level 
observables only, since the loop-induced observables are modified by new physics and are 
taken into account by the scanning procedure as discussed in the previous section. The 
uncertainties of the hadronic input parameters in the right-hand column are used to deter- 
mine the overall theory uncertainty on each observable, which then enters the x 2 function, 
added in quadrature with the experimental uncertainty. All uncertainties are assumed to 
be Gaussian. 

For the two most significant constraints, slightly looser uncertainties where assumed. 
For | |, the theory uncertainty was assumed to be at least 30% of the central value of 
the SM prediction. In the case of Rb, the partial decay width of the Z to b quarks, fits 
currently show a 2.5a deviation from the SM prediction [58], which is always worsened 
(or stays the same) in the models at hand [20]. In order not to be too restrictive, Rb was 
treated in the x 2 as if the experimental central value coincided with the SM central value, 
effectively requiring Rb to be within la of the SM value. Otherwise even SM-like points 
would be disfavoured. The same was done for the hadronic width R^, relevant in the U (2) 3 
case. 
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